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ABSTRACT 

A study  was  carried  out  to  determine  the  optimum 
placement  and  volume  of  a reinforcing  frame  around  a cut- 
out in  an  axially  loaded  stringer  and  ring  and  stringer 
stiffened  cylindrical  shell.  The  problem  was  analyzed 
using  the  linear  bifurcation  portion  of  STAGS  (Structural 
Analysis  of  General  Shells).  Four  parameters  were  varied; 
stringers  versus  rings  and  stringers,  cutout  size,  ratio 
of  frame  volume  to  cutout  volume,  and  frame  position.  It 
appeared  that  in  most  cases  the  position  with  the  frame 
next  to  the  cutout  edge  was  the  most  effective.  This  can 
be  attributed  to  the  frame's  ability  to  delay  the  onset  of 
local  buckling.  However,  there  was  a relative  maximum  in 
the  frame  distance  versus  critical  load  curves  for  a frame 
positioned  away  from  the  cutout  edge  at  a low  ratio  of 
frame  to  cutout  volume. 
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CUTOUT  REINFORCEMENT 


OF 

STIFFENED  CYLINDRICAL  SHELLS 
I . Introduction 

Background 

The  cylindrical  shell  is  an  important  aerospace  struc- 
tural element.  The  semi-monocoque  airplane  shell  structure 
of  the  1930' s was  typified  by  a very  thin  skin  supported 
by  much  stiffen  longitudinal  and  tranverse  reinforcing  mem- 
bers. As  a result,  buckling  of  the  thin  skin  panels  was 
primarily  a local  phenomenon.  With  the  advent  of  missies 
and  the  increase  in  flight  speed,  aircraft  skins  have  be- 
come thicker,  and  skins  and  reinforcement  are  of  similar 
rigidity.  Buckling  changed  from  a local  characteristic  to 
a problem  of  general  shell  instability.  A further  compli- 
cation is  encountered  when  an  access  door  or  similar  hole 
must  be  inserted  in  the  structure,  thus  weakening  an  other- 
wise very  efficient  structure.  The  purpose  of  this  thesis 
is  to  determine  the  strengthening  effects  of  a reinforcement 
around  a rectangular  cutout  in  a stiffened  cylindrical  shell. 

First,  a history  of  the. buckling  of  cylindrical  shells 
is  presented  to  put  this  problem  in  its  proper  context. 

Lorenz  (Ref  1)  was  the  first  to  derive  a buckling 

i 1 


J 

i 


formula  for  an  axially  compressed  cylindrical  shell.  He 
assumed  that  the  edges  were  simply  supported  and  that  both 
the  prebuckling  and  buckling  deformations  were  axisymmet- 
ric.  His  equations  were  slightly  inaccurate,  however. 
Timoshenko  (Ref  2)  published  the  correct  classical  buckling 
equation  in  1910.  The  classical  buckling  formula  is: 

S , = ( 3( 1-nu2) )~^  E ( h/R  ) = .6E(h/R)  (1) 

cl 

where  Scl  is  the  classical  buckling  stress. 

The  general  case  of  buckling  without  restriction  to 
axial  symmetry  of  the  buckling  deformations  was  first  dis- 
cussed by  Lorenz  in  1911,  assuming  that  displacements  in 
an  axial  direction  were  negligibly  small  and  the  boundary 
conditions  were: 

v=w=Mx=0  (2) 

The  coordinate  system  used  in  this  thesis  is  shown  in  Fig.  1. 

The  buckling  displacements  are  u , v , and  w in  the 
axial,  circumferential  and  radial  directions,  respectively. 
The  coordinate  directions  are  x,  e,  and  z. 

In  1914,  Timoshenko  (Ref  2)  determined  that  the  non- 
axisymmetric  buckling  stress  was  equal  to  .6E(h/R),  the 
same  value  obtained  for  the  classical  axisymmetric  case. 

For  non-axisymmetric  buckling,  theory  predicts  the  entire 
surface  will  be  covered  by  relatively  small  buckles.  This 
buckling  pattern  will  be  referred  to  as  either  the  chess- 
board pattern  or  the  diamond  pattern. 

There  are  two  basic  problems  with  the  buckling  patterns 

predicted  by  the  linear  classical  theory. 
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The  first  is  that 
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Figure  1.  Shell  Geometry  and  Sign  Convention 
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axially  symmetric  buckling  can  occur  only  when  there  is 
some  plasticity,  while  classical  buckling  theory  assumes 
perfect  elasticity.  The  second  problem  is  that  the  uniform 
pattern  of  buckles  in  chessboard  buckling,  is  observed  in 
experiments  as  only  one  or  two  rows  of  buckles  with  the 
remainder  of  the  shell  smooth. 

In  1932  Flugge  (Ref  1)  tested  a number  of  celluloid 
cylinders  for  a variety  of  length  to  radius  ratios  and 
radius  to  wall  thickness  ratios.  lie  discovered  that  these 
buckled  at  stresses  considerably  below  those  predicted  by 
the  classical  buckling  equation,  with  experimental  to  theo- 
retical buckling  stress  ratios  ranging  from  .52  to  .65. 

To  determine  the  possible  causes  of  the  discrepancy  between 
experiment  and  theory,  Flugge  examined  the  effect  of  bound- 
ary conditions  and  small  geometric  imperfections  upon  the 
buckling  stress. 

In  the  classical  approach  to  shell  instability,  the 
shell  was  assumed  to  expand  uniformly  inthe  radial  direction 
along  the  shell's  entire  length  due  to  the  applied  axial  load. 
Flugge  solved  for  the  axially  symmetric  deformations  of  a 
shell  whose  end  sections  are  restrained  in  the  radial  direc- 
tion. He  found  that  a curved  region  existed  along  the  merid- 
ian near  the  ends.  Increased  loading  caused  the  curved 
regions  to  extend  toward  the  middle  and  eventually  the  radial 
displacements  became  infinitely  large  at  the  critical  load. 

He  noted  that  the  curvature  caused  local  stress  to  exceed  the 

elastic  limit  even  if  S was  well  below  the  elastic  limit.  In  his 
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investigation  of  a sine  wave  deviation  from  perfect  circu- 
lar cylindrical  form,  Flugge  determined  that  the  imper- 
fection would  dominate  the  buckling  and  grow  infinitely 
large  as  the  classical  load  was  approached. 

Another  investigation  into  the  effect  of  boundary  con- 
ditions was  given  by  Hoff  in  (Ref  3)  and  outlined  in  (Ref  1). 
Starting  from  equations  by  Nachbar  for  a pressurized  spheri- 
cal  shell,  Hoff  was  able  to  determine  that  a semi-ir.f ini te 
shell  with  the  near  edge  perfectly  free  to  displace  would 
buckle  axisymmetrically  at  one-half  the  classical  buckling 
stress.  In  a follow  up  paper  by  Nachbar  and  Hoff,  scr/sci 
was  determined  to  be  .38  for  non-axisymmetric  buckling. 

For  many  years,  the  biggest  problem  with  the  solution 
of  shell  buckling  was  the  complexity  of  the  governing  shell 
equations.  In  1933,  Donnell  (Ref  4)  formulated  equilibrium 
equations  for  a cylindrical  shell  under  torsion,  based  on 
the  following  assumptions: 

1.  The  material  is  perfectly  elastic. 

2.  The  thickness  is  small  compared  to  the  radius. 

3.  The  deflections  are  small  compared  to  the  thickness. 

4.  Neglect  distortion  due  to  transverse  shear. 

A more  convenient  form  of  Donnell's  small  displacement 
equations  (Ref  5)  allowing  general  loading  are: 

Rtvlx,x  + Nx  e,e  = Px 

RNxe,x  + Ne>,e  = pe  ^ ^ 

D V 4w  + Ne/R  = pr 
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The  applied  loads  are  p^ , pQ,  and  pr  in  the  axial, 
circumferential,  and  radial  directions,  respectively. 

The  Donnell  equations  have  the  advantage  that  they  are 
much  simpler  than  most  others  and  yet,  they  are  sufficiently 
accurate  in  buckling  problems  if  the  number  of  waves  around 
the  circumference  is  sufficiently  large.  In  1S47,  Batdorf 
(Ref  6)  used  these  equations  to  solve  a variety  of  stability 
problems,  and  compared  the  solutions  with  other  theoretical 
solutions  and  test  results. 

Using  the  Donnell  equations,  Hoff  (Ref  1)  demonstrated 
that  the  appropriate  boundary  conditions  for  the  classical 
case  of  buckling  are: 

“*  ■ w’-x*  - - v’  - ‘ 

This  is  not  the  only  possible  simple  support  condition. 

The  four  simple  support  conditions  proposed  by  Hoff  are 
(Ref  1): 
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In  this  notation  classical  simple  supports  correspond  to  SS3. 
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In  the  next  few  years  researchers  investigated  both  serai - 
infinite  and  finite  length  shells  with  various  edge  condi- 
tions with  the  aid  of  Donnell's  small  displacement  equations; 
and  the  results  compared  well  with  other  accepted,  but  tedi- 
ous solutions. 

The  significance  of  these  results  indicated  that  there 
were  solutions  other  than  the  classical  one  for  the  buck- 
ling of  a thin  walled  circular  cylinder  under  axial  com- 
pression, despite  the  fact  that  the  problem  was  governed  by 
linear  equations.  Another  important  aspect  these  solutions 
pointed  out,  w as  that  a thin  shell  placed  between  the  platens 
of  a testing  machine  can  buckle  anywhere  between  one-half 
and  the  full  value  of  the  classical  buckling  stress.  The 
exact  value  depends  on  the  friction  between  the  test  arti- 
cle and  the  platen. 

The  difference  between  the  predictability  of  the  buck- 
ling stress  of  bars  and  plates  and  thin  shells  had  puzzled 
Von  Karmen.  He  observed  that  buckles  grew  gradually  in  bars 
and  plates,  but  snap  through  in  shells.  After  buckling, 
bars  and  plates  carried  the  load;  while  shells  could  only 
carry  a fraction  of  the  buckling  load.  He  determined  that 
the  key  lay  in  the  curvature  of  the  shell.  The  classical 
solution  implies  that  the  surface  may  buckle  inward  or  out- 
ward with  equal  ease.  However,  real  shells  show  a prefer- 
ence to  buckle  inward.  In  1934  Donnell  (Ref  7)  remedied 
this  problem  by  adding  terms  that  were  non-linear  in  the 
displacements  which  took  into 
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account  the  non-linear  effects  of  the  shell  curvature 
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Note  that  the  linear  Donnell  snail  displacement  Eqs. 

(3)  may  be  obtained  from  the  non-linear  Donnell  large  dis- 
placement Eqs.  (6)  by  the  simple  omission  of  all  terms  non- 
linear in  u,  v,  and  w. 

The  problem  with  these  equations  was  that  they  were 
non-linear  partial  differential  equations  and  could  not  be 
solved  exactly.  Both  Von  Karmen  and  Donnell  obtained  rig- 
orous solutions  for  these  equations  by  assuming  displacement 
patterns.  In  1942,  Tsien  suggested  replacing  the  buckling 
stress  with  the  minimum  equilibrium  stress  after  buckling 
as  being  reasonable  for  buckling  purposes. 

While  others  were  working  on  physical  and  mathemati- 
cal reasons  for  the  strange  behavior  of  axially  compressed 
cylindrical  shells,  Yoshimura  (Ref  1)  found  a geometric 
reason  for  this  behavior.  In  1951,  Yoshimura  determined 
that  a circular  shell  could  buckle  into  a diamond  pattern 
similar  to  the  previously  mentioned  chessboard  pattern 
without  causing  any  membrane  stresses.  The  significance 
of  this  pattern  is  that  a thin  shell's  extensional  stiff- 
ness is  much  greater  than  its  bending  stiffness.  The 
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incompatibility  of , the  diamond  pattern  with  the  support 
condition  probably  is  the  reason  why,  experimentally, 
shells  buckle  with  only  one  or  two  rows  of  diamonds. 

In  1963,  Almroth  (Ref  1)  using  the  computer  investi- 
gated combinations  of  the  various  terms  of  the  series: 

w = I I A.,  cos  ( j 7T  x /L  ) cos  (k7Te  /L  ) (7) 

j=0  k=0  x y 

In  this  investigation,  Almroth  minimized  the  total  poten- 
tial energy  with  respect  to  mu  and  N where: 

mu  = m/n  N = n^(h/R)  (8) 

m is  the  number  of  half  waves  in  the  axial  direction  and  n 

. is  the  number  of  full  waves  in  the  circumferential  direction. 

1 

Almroth' s results  were  in  excellent  agreement  with  work  by 
Thielemann  (Ref  1),  who  had  performed  experiments  with 
extreme  care  and  later  compared  the  experimental  results 
with  analytical  results. 

The  problem  of  considering  various  combinations  of 

terms  in  Eq.  (7)  was  looked  at  by  several  authors  between 

1941  and  1965.  The  minimum  S R/Eh  determined  from  these 

cr 

efforts  was  .0427,  which  was  well  below  the  classical  value 
of  . 6 . 

It  is  interesting  to  note  that  with  the  larger  number 
of  terms  retained  in  the  expression  for  w , the  coefficients 
of  the  double  Fourier  series  approach  those  for  a Fourier 
series  expansion  of  the  Yoshimura  buckling  pattern.  A prob- 
lem with  the  technique  whereby  the  potential  energy  was  mini- 
mized  with  respect  to  the  wave  length  parameters  was  that, 
as  the  number  of  terms  increased,  the 
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solution  was  only  valid  for  an  infinitely  thin  shell,  since 
N was  found  to  approach  zero.  This  was  overcome  by  requir- 
ing N^4n/R,  since  n,  the  number  of  circumferential  waves 


could  not  be  less  than  two. 

Investigations  based  on  the  Donnell  large  displace- 
ment equations,  were  carried  out  by  several  authors  in  the 
next  few  years  to  determine  analytically  the  effect  of 
small  deviations  from  the  nominal  geometry  upon  the  buck- 
ling load  of  a thin  cylindrical  shell.  It  was  noted  by 
Madsen  and  Hoff  (Ref  1)  that  a non-symmetric  deviation 
amounting  to  one-tenth  the  wall  thickness  coupled  with  an 
amplitude  of  axisymmetric  deviations  one-fortieth  the 
thickness  of  shell  resulted  in  a reduction  of  maximum  load 
to  60  per  cent  of  the  classical  value  for  a perfect  shell. 

Hoff  noted  that  the  large  displacement  Donnell  equa- 
tions were  completely  inadequate  to  represent  the  inexten- 
sional  deformations  of  Yoshimura  pattern  when  there  are 
five  to  ten  triangles  around  the  circumference  of  the  shell. 
Mayers  ft  Rehfield  (Ref  1)  also  noted  that  the  equations 
could  not  handle  the  plastic  stresses  along  the  triangle 
edges  with  the  result  that  the  curvature  predicted  along 
the  triangle  edges  was  too  large. 


Another  limitation  of  the  theory  is  that  it  does  not 
adequately  express  the  reduction  of  the  buckling  stress  due 
to  large  R/h.  Theory  predicts  that  k is  independent  of  R/h 


small  even  using  the  Donnell  large  displacement  equations. 

Madsen  & Hoff  (Ref  1)  demonstrated  that  the  replacement 
of  the  grossly  inadequate  large  displacement  Donnell  equa- 
tions by  more  accurate  equations  does  not  significantly 
change  the  equilibrium  stresses  obtainable  from  the  equations. 
The  explanation  of  the  apparent  paradox  is  that  the  large 
number  of  buckles  around  the  circumierence  allow  the  shell 
to  be  considered  a shallow  shell. 

Stein,  Fischer,  and  Almroth  separately  investigated  the 
effect  of  prebuckling  displacements  on  the  bulling  stress 
for  a variety  of  boundary  conditions.  Almroth  compared  the 
results  and  concluded  that  the  effect  of  boundary  conditions 
was  much  larger  than  the  effect  of  prebuckling  displacements. 

Koiter  was  the  first  to  call  attention  to  the  importance 
of  the  stability  of  the  system  at  the  bifurcation  itself. 

He  demonstrated  that  if  the  second  derivative  of  the  poten- 
tial energy  was  zero,  and  if  the  third  derivative  were  not, 
then  the  system  would  be.  unstable.  This  is  true  since  the 
bifurcation  point  would  correspond  to  a minimax  of  the  poten- 
tial energy.  Under  such  circumstances  the  system  is  very 
sensitive  to  small  initial  deviations  from  the  perfect  shape. 

Stein  (Ref  8)  presented  an  excellent  survey  of  advances 
in  shell  buckling  during  the  time  period  of  1960-1967.  In 
this  paper,  Stein  dealt  with  two  maior  topics:  The  impact 

of  improved  non-linear  shell  prebuckling  theory  and  the 
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significance  of  experiments  upon  near-perfect  shell  speci- 
mens. 

Stein  gave  a brief  derivation  of  consistent  shell  buck- 
ling equations  to  point  out  the  differences  between  con- 
sistent and  conventional  buckling  theory.  In  conventional 
theory, linear  membrane  theory  is  assumed  to  apply  in  the 
prebuckling  range  which  leads  to  constant  radial  displace- 
ments. This  in  turn  assumes  that  the  shell  boundary  con- 
ditions change  during  buckling.  On  the  other  hand,  in 
consistent  buckling  theory  the  non-linear  large  deflection 
Donnell  equations,  Eqs.  (6),  are  used  so  that  the  prebuck- 
ling boundary  conditions  can  be  consistent  with  the  boundary 
conditions  for  the  buckling  solution.  It  should  be  noted 
that  the  buckling  equations  for  consistent  theory  involve 
prebuckling  terms  which  would  vanish  in  conventional  theory. 

Stein  picked  results  from  several  researchers  who  had 
tested  near-perfect  shells.  The  materials  used  were  photo- 
plastic, nickel,  and  copper.  The  importance  of  these  tests 
was  that  the  near-perfect  shells  buckled  at  stresses  con- 
siderably closer  to  the  classical  buckling  stress  than 
previous  experiments.  The  conclusion  was  that  most  of  the 
scatter  in  experimental  results  was  due  to  small  initial 
deviations  from  perfect  shell  geometry. 

The  use  of  non-linear  bending  theory  has  lead  to  closer 
agreement  between  theory  and  experiment,  particularly  for  a 
near-perfect  shell  with  R/h  less  than  300;  see  figure  1 
(Ref  8). 

12 


Stein  also  includes  a section  on  shallow  spherical 
shells  under  external  pressure  with  similar  improved  agree- 
ment between  theory  and  experiment. 

Practical  shells  do  not  have  near-perfect  geometry  as 
assumed  by  theory  and  therefore,  techniques  to  determine 
the  reduction  in  buckling  strength  due  to  imperfections 
are  required.  Artificial  buckling  criteria  such  as  post- 
buckling  minimum  loads  and  loads  corresponding  to  equal 
energy  states  do  not  account  for  imperfections,  and  are 
thus  inadequate. 

A measure  of  a shell's  sensitivity  to  initial  imper- 
fections is  given  by  the  imperfection  sensitivity  angle  s 
(Ref  8).  This  parameter  is  related  to  the  initial  post- 
buckling  slope  of  the  load  vs  characteristic  displacement 
curve,  and  is  defined  in  Fig.  2 taken  from  (Ref  8).  The 
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parameter  s varies  from  -3  to  1,  that  is  from  a condition 
where  it  doubles  back  on  itself  to  a condition  where  it 
is  tangent  to  the  prebuckling  curve.  For  s positive,  im- 
perfections do  not  appreciably  alter  the  buckling  load; 
for  s negative,  imperfections  will  decrease  the  buckling 
load.  s also  provides  a measure  of  the  expected  violence 
of  the  buckling  process;  for  s<-2  "snap  through"  buckling 
will  always  occur,  whereas  for  -2<s<0  buckling  will  be 
violent  or  gentle  depending  on  whether  loading  or  charac- 
teristic displacement  is  controlled. 

For  the  cylinder  in  axial  compression  the  sensitivity 
angle  s (calculated  on  the  basis  of  conventional  theory) 
is  s = -3.  This  condition  agrees  with  snap  through  buck- 
ling observed  in  experiments. 

Although  it  has  been  previously  pointed  out  that  the 
major  cause  of  disagreement  between  experiment  and  theory 
is  the  presence  of  initial  defects,  research  indicates 
that  boundary  conditions,  eccentricity  of  stiffening  and 
eccentricity  of  load  application  also  affect  the  buckling 
load . 

The  eccentricity  of  stiffeners  also  strongly  affects 
the  buckling  load.  For  example,  Cara  (Ref  1)  demonstrated 
that  externally  stiffened  shells  have  carried  more  than 

twice  the  load  by  internally  stiffened  cylindrical  shells. 
This  effect  will  be  covered  more  extensively  in  a later 
section.  , 5 ■ 

7 ' * 

Another  factor  which  aftects  the  buckling  load,  is 
the  eccentricity  of  the  loading  relative  to  the  neutral 
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axis.  Results  obtained  by  Block  and  Almroth  (Ref  1)  show 
that  up  to  four  times  the  load  are  available  with  the  load 
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applied  on  the  inside  of  the  neutral  surface  for  both  inter- 
nal and  external  stiffening  compared  to  the  cases  for  load- 
ing at  the  neutral  surface. 

An  aspect  of  stiffened  cylindrical  shells  that  is  not 
present  in  the  unstiffened  shell,  is  that  small  meridianal 
curvature  drastically  affects  the  buckling  load.  For  in- 
ternal stiffening,  the  addition  of  7.8  per  cent  outward 
rise  results  in  an  increase  in  load  by  a factor  of  nine. 
Unstiffened  shells  in  compression  do  not  carry  appreciably 
more  load  with  the  addition  of  small  meridianal  curvature. 

The  conclusion  that  can  be  drawn  is  that  properly  posed 
buckling  theory  is  correct;  and  that  there  is  no  need  for 
postbuckling  minimums. 

In  the  application  of  thin  shell  technology  for  aero- 
space applications,  it  is  often  necessary  to  design  a 
cylindrical  shell  with  a cutout  for  non-structural  access 
doors  to  internal  equipment,  or  any  of  a number  of  other 
purposes.  The  exact  solution  to  such  a problem  is  beyond 
present  analytical  capabilities.  As  a result,  several 
researchers  have  undertaken  projects  in  this  area.  One  of 
these  researchers  is  Starnes  (Ref  9),  who  has  investigated 
the  buckling  of  a thin  unstiffened  cylindrical  shell  with 
a single  cutout,  both  experimentally  and  theoretically. 

Starnes  (Ref  9)  conducted  two  series  of  experiments, 
one  on  shells  made  of  Dupont’s  Mylar  with  a lap  seam;  and 
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the  other  series  of  tests  on  seamless  electrof ormed  copper 

shells.  The  parametric  ranges  were:  400<R/h<960  and 

0<a/R<.5;  where  a was  the  hole  radius.  The  experimental 

buckling  loads,  displacements,  and  the  stress  distribution 

applied  at  the  end  of  the  shell  were  correlated  with  a 

theoretical  parametric  study  performed  by  means  of  a 

Ray leigh-R i tz  type  of  approximation.  The  results  of  these 

experiments  led  to  the  conclusion  that  the  governing  para- 

2 

meter  of  the  problem  was  related  to  a /Rh  rather  than  a/R 
as  suggested  by  Tennyson.  Lekkerkerker  (Ref  9)  has  shown 
that  the  parameter: 

mu'  = h ( 12 ( 1-nu2  ) (a^Rh)*8  (9) 

governs  the  prebuckling  stress  distribution  and  displace- 
ments for  a circular  cylinder  with  a circular  hole  in  its 
side.  It  should  be  noted  that  the  effect  of  the  hole  is 
very  localized  in  nature.  Both  membrane  and  bending  stress 
increments  occur;  but  the  bending  stress  increment  is  always 
much  less  than  the  membrane  stress  increment.  It  should 
be  appreciated  that  the  maximum  membrane  stress  will  rise 
significantly  above  the  stress  value  obtained  for  a hole  in 
a flat  plate. 

The  buckling  analysis  referred  to  later  and  the  experi- 
mental results  presented  below  indicate  that  the  buckling 

load  of  a circular  cylinder  with  a circular  hole  in  its 

2 

side  are  related  to  the  ratio  a /Rh.  It  is  reasonable  to 
assume  that  mu  of  Eq.  (9)  is  the  governing  parameter  for 
the  problem. 


16 


It  was  possible  to  identify  approximate  ranges  of 
the  parameter  mu'  with  different  buckling  characteristics. 
For  mu'  ^ .4  imperfections  were  the  dominating  factor 
causing  buckling  below  the  classical  load.  For  values  of 
mu'  between  .4  and  1.0  the  buckling  loads  dropped  sharply 
as  mu'  increased.  In  this  range,  the  prebuckling  stress 
concentration  around  the  hole  caused  local  buckling,  and 
since  the  entire  shell  was  sensitive  to  small  disturbances 
in  this  load  range,  general  buckling  was  caused  by  the 
local  snap  through.  For  mu'  between  1.  and  2.0  tnere  is  a 
transition  between  the  sharp  decline  in  buckling  load  of 
the  previous  region  and  the  milder  decline  beyond  mu'  = 

2.0.  For  mu'  greater  than  2 local  buckling  occurred  at  a 
low  enough  load  so  that  the  shell  was  not  sufficiently  sen- 
sitive to  disturbances  to  cause  general  collapse.  The 
collapse  load  was  in  general  slightly  higher  than  the  local 
buckling  load  for  mu'  > 2.  For  this  case,  the  maximum  dis- 
placement at  the  edge  of  the  hole  was  on  the  order  of  .25 
inch,  which  is  many  times  the  wall  thickness,  and  sur- 
prisingly large  when  compared  to  the  prebuckling  displace- 
ments in  a cylinder  without  a cutout. 

Starnes  (Ref  9)  presents  results  from  his  copper  shell 
experiments  for  the  range  0 £ mu'  £ 3.  The  resulting  curve 
of  S / S c ^ vs  mu'  followed  very  nicely  with  the  data  from  the 
Mylar  shells.  Of  interest  are  figures  6 and  7 in  (Ref  9) 
representing  the  prebuckling  displacements  at  S/Sc^  = .393. 
Starnes  performed 
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a simplified  analysis  of  the  buckling  of  a circular  cylin- 
drical shell  with  a single  circular  hole  assuming  that 
general  collapse  was  initiated  by  local  buckling  and  that 
the  stress  distribution  in  a flat  plate  closely  resembled 
the  membrane  stress  distribution  in  a cylinder.  Despite 
these  restrictive  assumptions,  the  trend  indicated  by  the 
analysis  followed  the  experimental  results  for  mu'  <'2.5. 

For  mu'  >2.5  the  Donnell  approximations  of  linear  shell 
theory  were  no  longer  valid  due  to  the  large  prebuckling 
displacements  and  theory  and  experiment  diverged  sharply. 

As '.previously  mentioned,  local  buckling  is  rarely 
critical  in  modern  stiffened  shells  and  hence  general  in- 
stability is  the  dominant  criterion.  For  many  years  stab- 
ility analyses  of  stiffened  shells  considered  equivalent 
orthotropic  shells  and  the  eccentricity  of  the  stiffeners 
was  not  taken  into  account  in  these  analyses,  due  to  the 
complexity  of  the  equations.  More  recently  with  the  aid  of 
computers  the  effect  and  influence  of  stiffener  eccentricity 
has  been  explored  by  authors  for  internal  and  external 
stiffening  for  both  rings  and  stringers.  The  general  con- 
clusion that  can  be  drawn  is  that  for  an  axially  loaded 
cylinder,  external  stiffening  is  most  efficient.  A simpli- 
fied approach  to  the  problem  of  incorporating  eccentric 
stiffening  into  the  shell  equations  is  called  smeared 
stiffener  theory  in  which  the  stiffeners  are  assumed  to  be 
evenly  distributed  or  "smeared"  over  the  entire  shell  sur- 
face. This  approach  appears  to  be  satisfactory  for  closely 
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stiffened  shells  that  fail  in  general  instability,  since 
the  effect  of  the  discreteness  of  stiffeners  is  usually 
negligible.  The  adequacy  of  smeared  stiffener  theory,  is 


still  in  question  and  must  be  verified.  A review  of  ex- 
perimental work  to  assess  the  bounds  and  extent  of  appli- 
cability of  linear  smeared  theory  was  given  Singer  (Ref  10). 

In  outline,  the  analysis  employed  the  Donnell  stability 
equations  through  the  force  and  moment  expressions.  In  the 
mathematical  model,  the  stiffeners  are  smeared  to  form  a 
cut  layer  of  many  parallel  rings  or  stringers  which  touch 
each  other  but  are  not  connected.  The  main  assumptions  are: 

a.  The  stiffeners  are  distributed  over  the  whole 
surface  of  the  shell. 

b.  The  normal  strains  vary  linearly  in  the  stiffener 
as  well  as  in  the  skin.  The  normal  strains  in 
the  skin  and  the  stiffener  are  equal  at  their 
point  of  contact. 

c.  The  inplane  shear  membrane  force  is  carried 
entirely  by  the  shell. 

d.  The  torsional  rigidity  of  the  stiffener  is  added 
to  that  of  the  skin. 

The  analysis  is  outlined  in  (Ref  10)  and  given  in  detail  in 
(Ref  11).  The  resulting  equations  have  been  solved  for  a 
variety  of  boundary  conditions.  For  clamped  ends  a Galerkin 
technique  was  employed  since  the  equations  could  not  be 
solved  directly. 

The  main  objective  of  investigating  linear  smeared 
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stiffener  theory  is  to  determine  the  effectiveness  of  the 
stiffeners  and  their  relative  importance.  Batdorf  (Ref  6) 
has  shown  that  the  shell  geometry  can  be  characterized  by 
a single  parameter  Z. 

Z . (1-nu2)^  (L2/Rh)  (10) 

The  important  stiffener  parameters  are  spacing,  shape,  cross 
sectional  area,  and  eccentricity. 

Gross  sectional  area  is  usually  the  most  important 
geometric  property  of  the  stiffener  determining  their  effect- 
iveness. However,  it  also  directly  affects  the  weight  and 
must  be  minimized.  For  an  axially  loaded  cylindrical  shell, 
stringers  are  the  most  efficient  stiffeners.  Outside 
stringers  have  been  shown  to  be  always  more  efficient  than 
an  equivalent  weight  unstiffened  shell. 

Rings  are  much  less  effective  than  stringers  in  the 
stiffening  of  cylindrical  shells  under  axial  load,  when 
alone.  But,  when  they  are  added  to  an  axially  stiffened 
shell,  the  buckling  load  is  increased  significantly. 

As  was  done  for  unstiffened  cylindrical  shells,  the 
effects  of  boundary  conditions  have  been  investigated  for 
stiffened  cylindrical  shells.  The  effects  for  external  and 
internal  stringers  are  significantly  different,  and  depend 
strongly  on  the  stringer  geometry.  In  general,  the  axial 
restraint  (u*  = 0 instead  of  Nx  = 0)  is  the  predominant 
factor  for  simple  supports.  With  internal  stringers,  axial 
restraint  (u*  = 0 instead  of  Nx  = 0)  increases  the  buckling 
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load  about  45% 


With  external  stringers,  axial  restraint  is  less  effec- 
tive, and  for  weaker  stringers  the  ef f ecti veness  of  axial 
restraint  decreases  even  more. 

For  stringer  stiffened  shells  under  axial  compression, 
rotational  restraint  is  a major  factor  which  raises  the 
buckling  load  considerably. 

The  effect  of  non-linear  prebuckling  deformations  were 
considered.  Singer  (Ref  10)  found  that  non-1 inear  prebuck- 
ling deformations  apparently  were  not  a major  factor  in  the 
determination  of  the  buckling  load  of  stiffened  cylindrical 
shells.  However,  for  short  stringer  stiffened  cylindrical 
shells,  the  non-linearity  in  the  prebuckling  deformations 
may  warrant  further  study. 

Another  boundary  effect  of  stringer  stiffened  shells 
is  the  load  eccentricity.  The  usual  analyses  assume  that 
the  load  is  applied  at  the  skin  midsurface.  The  end  moments 
caused  by  this  eccentricity  may  have  a significant  effect 
on  the  buckling  load.  The  length  of  the  shell  affects  the 
magnitude  of  this  effect.  That  is,  for  a short  shell  this 
effect  is  small,  while  for  a longer  shell  this  effect  is 
much  greater.  Moments  which  tend  to  bend  the  shell  middle 
outward  (and  give  rise  to  tensile  hoop  stresses),  increase 
the  buckling  load.  Hence,  for  internal  stringers,  loading 
through  the  shell  midsurface  decreases  the  buckling  load. 

In  1962,  Professor  Van  der  Neut  stated  that  he  expected 
a reduction  in  imperfection  sensitivity  in  stringer  stiffened 
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shells  with  the  conclusion  that  linear  theory  would  be 
adequate  to  predict  their  buckling  loads  for  general  insta- 
bility. 

Singer  indicates  that  stiffener  area  ratio  is  an  impor- 
tant factor  in  stringer  stiffened  shells.  It  is,  however, 
not  the  only  important  parameter  in  the  problem;  the  shell 
geometry  parameter  2 is  of  equal  importance. 

Initial  imperfections  have  been  identified  as  the  main 
cause  of  the  large  discrepancies  between  experimental  and 
theoretical  buckling  loads  of  unstiffened  cylindrical  shells 
under  axial  compression.  For  closely  stiffened  shells,  the 
influence  of  imperfections  is  less  pronounced,  however,  im- 
perfections remain  the  primary  degrading  factor. 

Another  factor  recently  shown  to  cause  scatter  in  ex- 
perimental buckling  stress  is  the  effect  of  inelasticity  at 
low  values  of  R/h.  This  plasticity  occurs  even  when  Scr  is 
well  below  the  elastic  limit.  Hence,  variation  of  the 
elastic  limit  of  the  material  could  cause  some  scatter  at 
all  values  of  R/h. 

Singer  discusses  various  attempts  by  several  authors  to 
determine  the  optimum  structural  efficiency  of  stringer 
stiffened  cylindrical  shells  under  axial  compression.  It  is 
shown  that  the  most  efficient  stiffening  is  for  lighter 
stringers  and  thinner  shell  as  opposed  to  thicker  equiva- 
lent shell  or  heavy  stringers.  For  the  oj timui  fi 
.3<A.,/bh<  .8  for  a stringer  stiffened  shell,  where  A ..  is 
the  cross-sectional  area  of  the  stringer  and  b is  the 
stringer  spacing. 

Practically  all  minimum  'weight  design  studies  have  been 
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based  on  linear  theory.  However,  care  must  be  taken  in 
such  studies  because  the;-'*  are  cases  where  an  unstiffened 
shell  is  predicted  to  be  more  efficient  than  a slightly 
stiffened  shell.  Thus,  an  erroneous  conclusion  may  be  drawn 
from  such  a study  if  the  much  higher  imperfection  sensitivi- 
ty of  the  unstiffened  shell  is  not  taken  into  account,  since 
this  could  well  reverse  the  relation  of  the  actual  buckling 
loads. 

Palazottc  (Ref  12)  explored  the  validity  of  using 
linear  bifurcation  theory  for  the  buckling  analysis  of 
stringer  and  ring-stringer  stiffened  cylindrical  shells  with 
rectangular  cutouts.  This  analysis  was  performed  using  the 
STAGS  (Ref  13)  (Structural  Analysis  of  General  Shells)  com- 
puter program.  This  program  incorporates  linear  bifurca- 
tion theory,  non-linear  collapse  analysis,  discrete  stiffen- 
ing, and  smeared  stiffening  theory. 

The  buckling  loads  of  stringer  and  ring  and  stringor 
cylindrical  shells  were  computed  using  both  the  linear  bifur- 
cation analysis  and  the  non-linear  collapse  analysis  capa- 
bility of  STAGS.  When  the  results  for  both  the  linear  and 
non-linear  analyses  compared  within  a few  per  cent,  Palazotto 
concluded  that  linear  bifurcation  theory  adequately  predicts 
the  buckling  load  of  a stiffened  cylindrical  shell  with  cut- 
outs. The  importance  of  this  fact  is  that  a non-linear 
collapse  analysis  is  much  more  expensive  than  a linear  bifur- 
cation analysis  in  terms  ol  computer  time. 

A limited  investigation  was  carried  out  to  determine 
the  effect  of  cutout  reinforcement  of  stiffened  cylindrical 
shells. 
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It  is  the  purpose  of  this  thesis  to  determine  the  vari- 
ation of  buckling  load  with  the  location  and  volume  of  a 
reinforcing  frame  around  a rectangular  cutout  in  a stiffened 
cylindrical  shell. 

General  Procedure 

The  buckling  load  of  the  cylinders  investigated,  were 
determined  using  the  linear  bifurcation  capability  of  the 
STAGS  (Structural  Analysis  of  General  Shells)  computer 
program  and  smeared  stiffener  theory.  Four  parameters  were 
varied  in  this  study:  stringer  versus  ring  and  stringer 

shell  stiffening,  cutout  size,  reinforcing  frame  to  cutout 
volume  ratio,  and  frame  position. 


II.  Theory  and  Model inq 


General  Description  or  Staqs 


STAGS  (Structural  Analysis  of  General  Shells)  is  a 
computer  program  to  analyze  the  behavior  of  general  shells 
under  arbitrary  static  and  thermal  loading.  Non-linearities 
caused  by  plasticity  and  deformations  are  accounted  for. 
STAGS  is  based  upon  an  energy  formulation  in  which  the 
derivatives  are  replaced  by  their  two-dimensional  finite 
difference  approximations.  Minimization  of  the  potential 
energy  leads  to,  in  the  general  case,  a non-linear  set  of 
algebraic  equations  which  are  solved  by  a modified  Newton- 
Raphson  technique. 

STAGS  can  perform  either  elastic  or  plastic  analyses 
of  shells.  For  the  elastic  portion  of  STAGS  three  types 
of  analyses  may  be  performed;  linear  elastic,  classical 
bifurcation  with  linear  prebuckling  displacements,  and  non- 
linear elastic  collapse  analysis.  S'J'AGS  also  performs  a 
non-linear  inelastic  analysis,  but  does  not  permit  the 
temperature  or  material  properties  to  vary  with  the  space 
coordinates. 

Geometry  is  defined  using  any  one  of  the  standard  shell 

geometry  options  (plate,  circular,  cylinder,  etc.);  or  the 
user  can  easily  'write  a subroutine  to  define  his  own  partic- 
ular oeonetry. 

Shell  stiffening  can  be  accomplished  in  several  ways. 
First,  discretely  located  stiffeners  can  be  positioned 


J 


anywhere  on  the  shell  surface  along  a coordinate  line. 

Second,  the  entire  shell  may  be  stiffened  by  using  the 
"smeared"  stiffener  capability.  Lastly,  a subroutine  can 
be  added  to  the  program  to  define  the  stiffness  properties 
of  the  shell  wall,  to  which  stiffener  effects  can  be  added. 
The  user  also  has  the  option  to  specify  one  of  the  standard 
wall  configurations.  Provisions  have  been  made  to  incor- 
porate shell  cutouts  along  coordinate  lines. 

Outline  of  STAGS  Theory 

A review  of  the  theory  used  in  STAGS  is  given  below. 

A more  complete  treatment  of  the  theory  is  given  in  (Ref  14). 
i 

The  solution  procedure  used  in  STAGS  is  based  upon  the 
principle  of  stationary  potential  energy.  The  shell  sur- 
face is  covered  with  mesh  lines  parallel  to  the  coordinate 
lines.  The  degrees  of  freedom  are  the  displacements  u,  v, 
and  w in  the  axial,  circumferential,  and  radial  directions, 
respectively. 

The  total  strain  energy  of  the  shell  is  obtained  by 
numerically  integrating  the  strain  energy  density  over  the 
surface  of  the  shell.  For  the  integration  procedure  the 
shell  surface  is  divided  into  a number  of  subareas,  in 
which  each  subarea  corresponds  to  the  intersection  of  two 
mesh  lines.  The  displacements  and  their  derivatives  are 
replaced  by  the  appropriate  finite  airference  approximations. 
A vector,  Z1 , of  strains  and  curvature  changes  is  evaluated 
at  mesh  station  i,  where: 

| 

i 

i 

L 
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Z indicates  the  transpose  of  Z . Note  that  Z is,  in 
general,  a non-linear  function  of  the  displacements  and 
their  derivatives  at  mesh  station  i. 

The  displacements  and  derivatives  at  mesh  station  i 
are  approximated  using  finite  difference  techniques,  hence, 
the  approximated  displacements,  their  derivative,  and  Z1 
are  functions  of  the  "raw"  displacement  unknowns  surround- 
ing mesh  station  i.  Excellent  descriptions  of  finite 
difference  approximations  can  be  found  in  the  works  by 
Deschler  (Ref  16)  and  Wiley  (Ref  17),  therefore  finite 
difference  approximations  will  not  be  discussed  further. 

The  total  vector  of  raw  displacement  unknowns  is: 

'r  I 1 i i i 1max| 

where  u1,  v1 , and  w1  are  the  "raw"  displacement  unknowns  at 

mesh  station  i,  and  i is  the  total  number  of  mesh  stations 

max 

Note  that  there  are  3*i  displacement  unknowns. 

max  K 

The  strain  energy  density  at  mesh  station  i/^U1,  is  a 
function  of  the  strains  and  curvature  changes.  Let  the 
matrix  D1  at  station  i be  defined  as: 


( 13) 


where  a subscript  following  a comma  indicates  differenti- 
ation with  respect  to  one  of  the  strains  or  changes  of  curva- 
ture. For  k = 1,2, 3, 4, 6, 6,  the  derivatives  are  to  be  taken 


with  respect  to  e , e , e. 


K , K , respectivelv 

^ 7 v A 7 1 J 


A derivation  of  D1  for  a shell  with  smeared  stringers  is 
given  in  (Ref  13).  The  strain  energy  density  at  subarea 


ax  is : 


. T 

AU1  = Z1  DX  Z1 


(14) 


A li*  is  a fourth  degree  polynomial  in  the  displacement  un- 
knowns. Summing  strain  energy  contributions  for  all  sub- 
areas  yields  the  total  strain  energy  for  the  shell. 


U a I A U1  a1 


(15) 


If  F is  the  vector  of  external  forces  corresponding  to 
the  displacement  unknowns,  then  the  potential  energy  of  the 
work  done  by  external  forces,  VI,  can  be  expressed  as 


W = XTF 


(16) 


The  total  potential  energy  of  the  system  is: 


V = U-W 


(17) 


Note  that  the  potential  energy,  V,  is  a fourth  degree  poly- 
nomial in  the  displacement  unknowns.  Static  equilibrium 
requires  that  the  first  variation  of  the  potential  energy 
be  ecual  to  zero  and  leads  to  the  equation: 


LX  = F 


( ia ) 


•where  the  non-linear  "stiffness"  operator  for  the  displace- 
ment unknowns  is  defined  by  the  following  equation: 

= grad  U 

■jc 


LX 


(19) 


IL.XI  . , of  LX  are  defined  as: 
i 


or  equivalently  the  rows, 


[LXj  . = ,j  U ( 20) 

1 2 X( i 5 

t h 

where  X(i)  is  the  l component  of  trie  vector  of  displace- 
ment unknowns. 

When  only  linear  terms  are  included  in  the  strain  and 
curvature  chance  relations,  L degenerates  into  a linear 
operator  of  matrix  form. 

When  L is  non-linear,  iterative  techniques  must  be  em- 
ployed for  the  solution  of  Eq.  (13).  A general  collapse 
analysis  requires  that  Eq.  (18)  be  solved  for  a sequence  of 
applied  loads.  A reasonably  good  initial  approximation  of 
the  solution  and  a series  of  moderately  increasing  load 
steps  are  required  for  reliable  detection  of  collapse  aue 
to  the  non-uniqueness  of  the  solutions  to  non-linear  equa- 
tions. The  requirement  that  Eq.  (18)  be  solved  repeated) y 
using  iterative  techniques,  is  what  makes  a non-linear 
collapse  analysis  so  expensive  compared  to  a linear  bifur- 
cation analysis  (described  later). 

A brief  description  will  be  given  of  a modified  Newton 
Raphson  method  used  in  the  solution  of  Eq.  (18)  for  the  dis 
placement  unknowns,  X.  For  a function  g(x)  = 0,  the  well- 
known  Newton-F aphson  iteration  is: 


f 


(21) 


where  X.  is  an  approximate  solution  to  q(x)  = 0 and  X.  . is 
i l + l 

a better  approximation.  Repeated  application  of  Eq.  (21) 
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yields  a solution  to  g(x)  = 0 to  any  desired  accuracy 
Notation  will  be  simplified  with  introduction  of: 


2 U 


(22) 


1 3 2 xTi  )<)  X ( j ) 

» 

Unless  the  non-linear  terms  of  L are  dropped,  L will  be  a 
function  of  a particular  displacement  vector  X,  , and  will 


be  designated  L 


With  L 


the  Newton-Raphson  method 


k 

can  be  readily  extended  to  solve  Eq.  (18).  The  iteration  is: 


xk.i  " **  ♦ \ ' <F-L>:k> 


I _ -1  t 

Lv  * indicates  the  inverse  of  L,.  , 
Xk  Xk 


(23) 

The  inversion  requires 


that  Ly  be  non-singular,  that  is,  a unique  solution  of 


Xh 

v » 

Eq.  (13)  requires  that  the  determinant  of  L be  non-zero. 

k 


The  most  effective  use  of  Eq.  (23)  in  the  solution  of 


Eq 

of 

f 


(18)  suggests  only  periodic,  recomputation  and  ' inversion 

» 

Ly  , since  the  computer  execution  time  required  to  invert 
Ak 


Ly  is  high  compared  to  the  time  required  to  evaluate  the 

k 

reaminder  of  Eq.  (23).  The  modified  form  of  the  Newton- 
Raphson  method  used  in  STAGS  is: 


xk*i  • vC  <f-lv  <24) 

m 

where  the  subscript  m indicates  the  last  iteration  for  wnich 

» 

Lv  ' was  recomputed, 
m 

The  mathematical  characterization  of  bifurcation  buckling 
is  provided  by  the  generalized  Newton-Raphson  method.  Let 
X satisfy  Eq.  (18)  under  the  load  vector  F,  then  if  for 


\A 


■ 

' 


. 

■ 

I ; 


every  neighborhood  of  X there  exists  another  vector  Y such 
that 

LXY  = F (25) 


then  bifurcation  buckling  has  occurred  at  load  F.  As  pre- 
viously mentioned,  a unique  solution  of  Eq.  (18)  using  the 

I 

iteration  scheme  of  Eq.  (23)  requires  that  L v be  non- 

Xk 

singular.  Conversely,  if  there  are  multiple  solutions  of 

9 

Eq.  (18)  then  L v is  singular,  in  which  case: 

Xk 

det  (l'y  ) = 0 (26) 

Xk 

Let: 

X.XXL  (27) 


where  X is  the  vector  of  displacement  unknowns  that  satis- 
L/ 

fies  Eq.  (18)  for  a given  load  vector  F^.  Substituting  Eq. 
(27)  into  equation  (26)  we  obtain: 


det  (L 


0 


(28) 


For  classical  buckling  in  which  prebuckling  rotations  are 

f 

ignored,  only  the  linear  terms  of  L ^ x are  retained,  result- 

L 

ing  in  an  eigenvalue  problem  of  the  form: 


det  (A-AB)  = 0 


(29) 


The  specific  method  employed  to  obtain  th  "undamental 
eigenvalue  A and  its  eigenvector  is  not  mentioned  in  the 
STAGS  manuals  (Ref  13)  and  (Ref  14).  If  the  reader  wishes 
more  information  in  this  area,  a description  of  the  Power 
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A 


Method  can  be  found  in  (Ref  18)  and  several  methods  as 
applied  to  a large  finite-element  program  are  described  in 
the  NASTRAN  Theoretical  Manual  (Ref  19).  The  buckling 
load  vector  Fcr  is  a multiple  of  the  applied  load  vector 
for  the  linear  solution. 


F 


cr 


=*fl 


(30) 


The  formation  of  the  A and  B matrices  in  Eq.  (29) 
will  be  considered  briefly.  From  Eq.  (15)  we  know  that: 


b 2U  = T ak  32*uk  (31) 

jwnrnvf  k 3 x(  i )3x(  j ) 

Using  Eq.  (14)  the  kt*1  term  of  Eq.  (31)  is: 

32Uk  = 3Zk  Dk  3Zk  +3  2Zk  ASk  (32) 

3xd  )3X(  j ) JxTn  JxTjT  TxTTTTxTTT 

where  ST  is  the  linear  stress  resultant  vector  at  mesh 
station  k.  The  more  general  (non-linear)  stress  resultant 
vector  is: 

Sk  = Dk  Zk  (33) 


where 

kT 

S = |N  N N M M M I (34) 

[_  x e xe  x e xej 

The  terms  generated  by  the  first  term  on  the  right  side 
of  Eq.  (32)  make  up  the  A matrix,  and  represent  contribu- 
tions from  membrane  energy. 

The  second  term  in  Eq.  (32)  contributes  to  the  B matrix 
and  represents  the  bending  energy. 
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An  example  of  the  derivation  of  the  constitutive  rela- 


tions, Eq.  (33),  for  a shell  with  smeared  stringers  is 
given  in  (Ref  13). 

Modeling 

The  shells  investigated  in  this  thesis  were  stiffened 
circular  cylindrical  shells  with  reinforcements  around  two 
symmetrically  placed  rectangular  cutouts.  The  linear  bifur- 
cation capability  of  STAGS  was  used  to  predict  buckling 
loads,  as  opposed  to  a non-linear  collapse  analysis  which 
would  have  required  considerably  more  computer  time.  The 
accuracy  of  linear  bifurcation  theory  for  stiffened  cylin- 
drical shells  using  STAGS  was  substantiated  by  Palazotto 
in  (Ref  12),-  He  analyzed  a stringer  stiffened  shell  with 
two  rectangular  cutouts.  The  finite  difference  mesh  con- 
sisted of  602  mesh  points,  with  a concentration  of  mesh 
points  in  the  vicinity  of  the  cutouts.  Due  to  triple 
symmetry  only  one-eighth  of  the  shell  was  analyzed. 

The  buckling  load  for  this  shell  was  computed  using 
both  linear  bifurcation  theory  and  non-linear  collapse 
theory.  Despite  the  fact  that  the  radial  displacements 
were  considerably  greater  for  the  non-linear  analysis,  the 
buckling  load  predicted  by  linear  bifurcation  theory  was 
only  6.87  per  cent  higher  than  the  collapse  load  predicted 
by  non-linear  theory.  This  is  an  acceptable  level  of  accur- 
acy since  the  expense  of  non-linear  analyses  would  severely 
limit  the  affordable  range  of  parameters  that  could  be 
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investigated. 

Smeared  stiffener  theory  was  used  to  represent  the  uni- 
formly spaced  rings  and  stringers  stiffening  the  entire 
shell.  An  alternate  method  was  to  specify  a discrete  stif- 
fener for  ring  or  stringer.  The  additional  input  time  re- 
quired for  discrete  stiffening  was  not  justified  by  a 
commensurate  increase  in  accuracy  of  the  buckling  load  as 
shown  by  Palazotto  (Ref  12).  He  indicated  that  the  dif- 
ference in  buckling  loads  using  smeared  stiffeners  and 
discrete  stiffeners  amounted  to  less  than  four  percent 
for  a similar  problem. 

On  the  basis  of  the  number  of  mesh  points  used  by 
Palazotto  and  current  geometry,  the  author  used  a minimum 
of  776  mesh  points. 

The  reinforcing  frame  around  the  cutout  was  modeled 
using  discretely  located  stiffeners. 
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III.  Results 


Model  Description 

In  this  thesis,  a stiffened  circular  cylindrical  shell 
with  a reinforcing  frame  around  the  cutout  was  analyzed  for 
a variety  of  parameters.  The  basic  geometry  was  held  con- 
stant while  four  parameters  were  varied  independently: 


Internal  Stiffening 

Stringers  Only 
Rings  and  Stringers 


Total  Cutout  Size  2a 

12  x 12  in.  total  cutout  size 
24  x 24  in.  total  cutout  size 

Ratio  of  reinforcing  frame  volume  to  removed  cutout 


volume 

W ,/W  = 

sf  o 

IT) 

• 

O 

W c/^i  = 
sf  o 

1.0 

w ,/w  = 

sf  o 

1.5 

W , /w  = 
sf  o 

2.0 

Reinforcing  frame  position  (4  positions) 

The  basic  shell  geometry  is  shown  in  Fig.  3.  The 
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Cutout  Sizes 
2a  = 12" 

2a  = 24" 


E = 10'psi 
nu  •=  .333 
G = 3.75* 106psi 
Z = 1591. 


Y 

' x 

Portion  of  Shell 
Analyzed 


Figure  3 


Easic  Shell  Geometry 


— — ”1 

radius  of  the  shell  was  chosen  such  that  the  circumferen- 
tial coordinate  was  numerically  equal  in  either  degrees  or 
inches.  The  length  was  chosen  such  that  a finite  difference 
mesh  could  be  imposed  upon  the  shell  compatible  with  the 
stiffener  and  reinforcing  frame  locations  and  the  cutout 
boundaries,  since  STAGS  requires  that  the  discrete  stiffeners 
and  the  cutout  edges  lie  along  mesh  lines.  Smeared  stiffe- 
ners should  also  lie  along  mesh  lines  for  accurate  repre- 
sentation. The  value  of  Z for  this  shell  is  1591.  The 
overall  dimensions  were  chosen  so  that  this  thesis  could  be 
compared  with  previous  work  by  Palazotto  (Ref  12).  The 
dimensions  given,  approximate  a large  missle  interstage. 

The  smeared  stiffener  geometry  and  spacing  are  indica- 
ted in  Fig.  4.  The  ring  and  stringer  geometries  were  de- 
cided upon  in  order  to  approximate  those  that  would  be 
expected  in  actual  applications.  The  two  to  one  ratio  of 
the  height  to  width  was  chosen  as  compromise  of  good  bend- 
ing stiffness  and  yet  reasonable  torsional  stiffness. 

The  two  cutout  sizes  used  were  a = 6"  and  a = 12" 

corresponding  to  12"  x 12"  and  24"  x 24"  cutouts  respect- 

u 

ively.  The  value  of  the  cutout  size  parameter,  a/(Rh)  , 
is  2.57  and  5.13  for  a = 6"  and  a = 12",  respectively. 

The  value  cf  a = 12  approximated  the  largest  cutout  size 
for  which  linear  theory  is  valid  (Ref  12). 

The  reinforcing  frame  volume  to  removed  cutout  volume 
ratios  correspond  to  the  situation  where  the  reinforcing 
frame  weighs  from  one-half  to  twice  the  weight  of  the 


> 


material  removed  for  the  cutout.  This  was  considered  as 

large  a range  as  would  be  of  practical  interest. 

Note  that  the  frame  positions  shown  in  Fig.  5 are 

different  for  a = 6"  and  a = 12".  The  maximum  frame 

positions  were  chosen  as  being  the  maximum  distance  at 

which  the  frame  would  have  a significant  impact  on  the 

buckling  load  as  determined  in  (Ref  12). 

Note  that  for  a given  volume  ratio  the  total  frame 

volume  remained  constant  as  the  frame  position  was  varied. 

This  was  necessary  in  order  to  validly  compare  the  most 

efficient  frame  position  for  a given  volume  ratio. 

The  equations  for  the  geometric  properties  of  the 

reinforcing  frame  will  be  discussed.  The  volume  per  unit 

area  of  the  cutout  for  the  stringer  stiffened  shell  is 

WOA  = .19068  IN^/IN^;  and  for  the  ring  and  stringer 

3 2 

stiffened  shell  WOA  =.20975  IN  /IN  . The  cross  sectional 
area  of  stiffening  frame  is  given  by: 


A f = a WOA  (Wsf/WQ)/2d 


(35) 


where  2a  is  the  total  cutout  dimension  and  d is  the  distance 
from  the  cutout  centerline  to  the  stiffening  frame.  For  a 
square  cross-section,  the  height  and  width  of  the  stiffen- 
ing frame  become: 


bsf  ■ hsf  - (Asf>  ' 


(36) 


The  moments  of  inertia  equal: 


bsf\f - hsf4'12 


(37) 
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The  torsional  stiffness  constant  (Ref  15),  for  b/h  = 1 
is 

J _ = .1406  h - (38) 

sf  sf  sf 

• For  the  present  value  of  G 

GJsf  = 527382.  hgf4  psi  (39) 

The  values  of  the  frame  properties  described  above  are 
given  in  Appendix  A.  The  listing  is  ordered  according  to 
computer  run  number  with  the  geometric  parameters  shown 
for  each  run. 

The  analysis  was  carried  out  over  one-eighth  of  the 
shell  surface  since  the  loading  and  shell  were  triply 
symmetric  about  the  planes  x = 48",  e = 0°,  and  z = 0". 
Sufficient  accuracy  with  the  minimum  number  of  mesh  points, 
required  separate  finite  difference  mesh  arrangements  for 
a = 6"  and  a = 12".  These  meshes  are  shown  in  Fig.  6. 

A uniform  axial  displacement  was  applied  at  the  end 
of  the  cylinder.  This  type  of  loading  used  is  comparable 
to  usual  experimental  load  application  (Ref  9).  Symmetric 
boundary  conditions  were  dictated  by  analyzing  the  problem 
over  one-eighth  of  the  shell.  The  boundary  conditions  and 
enforced  displacements  are  shown  in  Fig.  7. 
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iaure  6-  Finite  Di 


Numerical  Resul ts 


As  previously  stated,  a stiffened  cylindrical  shell 
with  a reinforcing  frame  around  a cutout  was  analyzed  for 
a variety  of  parameters  using  the  linear  bifurcation  por- 
tion of  STAGS. 


A tabular  listing  of  the  buckling  loads,  Pcr?  is  given 
in  Appendix  B,  according  to  computer  run  number.  In  addit- 
ion, the  ratios  P /P  and  P /P„  are  listed,  where  P and 

cr  o cr  1 o 

P^  are  dependent  upon  the  particular  shell  in  question. 

Po  is  the  buckling  load  of  a shell  with  the  same  type  of 

overall  shell  stiffening  (stringers  only  or  rings  and 

stringers)  as  the  shell  in  question,  but  one  that  does  not 

have  cutouts  or  reinforcing  frames.  P^  is  the  buckling 

load  of  a shell  as  described  above,  but  with  a cutout. 

The  buckling  load  ratio  P /P  is  plotted  versus  the 
y cr  o ^ 

distance,  d,  of  the  frame  from  the  centerline  of  the  cutout 
in  Figures  8 through  11.  Figures  8 and  9 are  for  stringer 
stiffened  shells;  and  Figures  10  and  11  are  for  shells  stif- 
fened by  rings  and  stringers.  For  Figures  8 and  10  the  cut- 
out size  is  2a  = 12";  and  for  Figures  9 and  10,  2a  = 24". 

P,,/P  has  been  shown  for  reference. 

1 o 

An  inspection  of  Figures  8 through  11  indicates  that 


the  most  effective  position  for  the  reinforcing  frame  is 
always  along  the  cutout  edge. 

It  has  been  observed  in  (Kef  12)  that  the  general 
collapse  shape  in  the  vicinity  of  a cutout  is  similar  to 
that  of  a shell’s  displacement  field  under  a linear  prebuck- 
ling analysis.  Thus,  it  is  possible  to  appreciate  the  ideal 
frame  position  by  comparing  the  linear  displacement  fields 


of  shells  with 
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the  reinforcing  frames  in  various  locations.  The  linear 
displacements  of  a stringer  stiffened  shell  with  a cutout 
size  2a  = 24"  and  volume  ratio  W ,/W  = 2,  are  shown  in 

Figures  12  and  13  for  the  case  with  no  frame  and  the  case 
with  a frame  at  the  cutout  edge,  respectively.  In  Figure 
12,  notice  the  large  radial  displacements  along  the  upper 
and  lower  edges  of  the  cutout.  In  Figure  13  the  very 
definite  stiffening  effect  of  the  frame  along  the  cutout 
edge  is  readily  apparent.  Both  the  large  radial  displace- 
ments along  horizontal  edges  and  the  large  rotations  along 
the  vertical  edges  of  the  cutout  have  been  dramatically 
reduced  by  the  reinforcement. 

An  additional  optimum  position  occurs  for  a stringer 
stiffened  shell  with  cutout  size  2a  = 12"  and  volume  ratio 
Wg^/Wo  = .5,  in  which  the  reinforcing  frame  is  equally 
effective  6"  away  from  the  cutout  as  it  is  along  the  cut- 
out edge.  This  is  in  spite  of  the  fact  that  at  the  outer 
position  the  cross-sectional  area  of  the  reinforcing  frame 
is  only  half  the  area  of  inner  frame  (frame  volume  is  inde- 
pendent of  position),  and  the  moments  of  inertia  have  a 
ratio  of  1:4.  The  linear  displacement  fields  are  shown  in 
Figures  14,  15  and  16  for  no  frame,  a frame  along  the  cutout 
edge,  and  a frame  6"  from  the  cutout  edge,  respectively. 
Inspection  of  Figures  14  and  15  reveals  that  the  displace- 
ments and  rotations  along  the  cutout  edge  are  again  some- 
what smaller  for  the  frame  adjacent  to  the  cutout  edge  than 
those  for  the  unreinforced  shell.  For  the  frame  6" 


away 


Shell  Orientation  in  the 
Displacement  Plots 


Linear  Displacements  For  a Stringer  Stiffened 
Shell,  2a  = 24" 


j 

: i 


RUN  29  LINEAR  X 25. 


Note:  The  frame  position  is  indicated  in  this  and 
subsequent  diagrams  by  a heavy  black  line 


Fiaure  13.  Linear  Displacements  For  Run  29 
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Linear  Displacements  For  a Stringer  Stiffened 
Shell,  2a  = 12" 
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Figure  15.  Linear  Displacements  For  Run  1 
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from  the  cutout,  however,  the  displacements  are  not  re- 
duced as  much.  Figures  17  through  19  display  the  axial 
compressive  stress  as  a surface  for  the  shells  mentioned 
above.  It  is  evident  that  a reinforcing  frame  in  posi- 
tion 1 significantly  reduces  the  stresses  in  the  vicinity 
of  a cutout  below  those  for  an  unreinforced  shell.  This 
stress  reduction  lessens  the  chance  for  local  buckling, 
consequently,  leading  to  a higher  shell  buckling  load. 
However,  even  a close  examination  of  Figures  17  and  19, 

.nows  that  a frame  at  position  4 yields  only  a slight  re- 
duction of  stress  in  the  area  of  the  cutout. 

Next,  a comparison  will  be  made  of  the  normalized 
displacements  (maximum  displacement  = 1)  occurring  at  bi- 
furcation; these  are  given  in  Figures  20  through  22.  Com- 
paring these  Figures  shows  that  the  buckling  displacements 
near  the  cutout  are  quite  different  for  the  reinforced  and 
unreinforced  shells.  Apparently,  even  the  small  frame  at 
position  4 is  sufficient  to  alter  the  buckling  mode. 

The  trends  of  the  curves  in  Figures  8 through  11  indi- 
cate that  a reinforcing  frame  positioned  far  away  from  the 
cutout  may  cause  the  shell  to  buckle  at  a load  lower  than 
the  buckling  load  for  an  unreinforced  shell.  This,  in  fact, 
occurs  in  Figure  10  for  a ring  and  stringer  stiffened  shell 

for  a volume  ratio  W £/W  = .5  v/ith  frame  at  position  6. 

sr  o 

In  order  to  suggest  an  explanation  for  this  occurrence, 
consider  a stringer  stiffened  shell,  with  2a  - 12", 

/WQ  = .5,  and  a frame  at  position  6.  The  plots  of  linear 
prebuckling  deflections  and  stresses  are  nearly  identical 
to  those  for  a shell  without  the  reinforcing  frame. 
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from  the  cutout,  however,  the  displacements  are  not  re- 
duced as  much.  Figures  17  through  19  display  the  axial 
compressive  stress  as  a surface  for  the  shells  mentioned 
above.  It  is  evident  that  a reinforcing  frame  in  posi- 
tion 1 significantly  reduces  the  stresses  in  the  vicinity 
of  a cutout  below  those  for  an  unreinforced  shell.  This 
stress  reduction  lessens  the  chance  for  local  buckling, 
consequently,  leading  to  a higher  shell  buckling  load. 
However,  even  a close  examination  of  Figures  17  and  19, 
shows  that  a frame  at  position  4 yields  only  a slight  re- 
duction of  stress  in  the  area  of  the  cutout. 

Next,  a comparison  will  be  made  of  the  normalized 
displacements  occurring  at  bifurcation;  these  are  given 
in  Figures  20  through  22.  Comparing  these  Figures  shows 
that  the  buckling  displacements  near  the  cutout  are  quite 
different  for  the  reinforced  and  unreinforced  shells. 
Apparently,  even  the  small  frame  at  position  4 is  suffi- 
cient to  alter  the  buckling  mode. 

The  trends  of  the  curves  in  Figures  8 through  11  indi- 
cate that  a reinforcing  frame  positioned  far  away  from  the 
cutout  may  cause  the  shell  to  buckle  at  a load  lower  than 
the  buckling  load  for  an  unreinforced  shell.  This,  in  fact, 
occurs  in  Figure  9 for  a ring  and  stringer  sitffened  shell 

for  a volume  ratio  W ,/W  = .5  with  frame  at  position  6. 

sf  o 

In  order  to  suggest  an  explanation  for  this  occurrance, 
consider  a stringer  stiffened  shell,  with  2a  = 12", 

Wgf/W  = *5,  and  a frame  at  position  6.  The  plots  of  linear 
prebuckling  deflections  and  stresses  are  nearly  identical 
to  those  for  a shell  without  the  reinforcing  frame. 
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Normalized  Buckling  Displacements  For  a 
Stringer  Stiffened  Shell,  2a  = 12" 
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Figure  21.  Normalized  Buckling  Displacements  For  Run  1 
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The  normalized  buckling  displacements  are  shown  in  Figure 


23.  Notice  that  for  this  frame  position,  the  buckling  pat- 
tern no  longer  resembles  the  patterns  for  frame  positions 
1 and  4;  instead,  the  pattern  is  similar  to  that  of  the 
shell  without  a frame.  Also,  notice  that  the  vertical  leg 
of  the  frame  lies  at  the  bottom  of  an  inward  buckle.  It  is 
supposed  that  the  eccentricity  effect  of  the  internal  frame 
contributes  to  the  early  buckling  of  the  shell  due  to  the 
additional  compressive  stress  caused  by  the  frame's  eccen- 
tricity leading  to  a less  stable  configuration. 

In  Figures  24  through  27  buckling  load  ratio  Pcr/p0  is 

plotted  versus  volume  ratio  W ,,/W  . The  variation  of  shell 

sf  o 

stiffening  and  cutout  size  is  in  the  same  order  as  for  Fig- 
ures 8 through  11.  As  a general  rule  the  heavier  reinforc- 
ing frames  (W  ,,/W  ) lead  to  higher  buckling  ratios  (P  /P  ). 
An  exception  to  this  trend  is  shown  in  Figure  26  for  frame 
position  6.  Notice  that  for  stringer  stiffened  shells, 
(Figures  24  and  25)  the  light  frames  (Wsf/WQ  = *5)  produce 
a strengthening  effect  that  is  roughly  independent  of  frame 
position.  The  effect  of  frame  position  does  not  become 
large  until  heavier  frames  are  used.  On  the  other  hand, 
for  ring  and  stringer  stiffened  shells  (Figures  26  and  27) 
the  effect  of  frame  position  is  very  large  for  light  frames 

( W f/W  = .5). 
s f o 

Apparently  the  light  frames  for  ring  and  stringer 
stiffened  shells  are  affecting  a factor  of  buckling  which 
; s not  significant  in  stringer  stiffened  shells.  An 
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P /P  vs  W For  a Stringer  Stiffened  Shell, 
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interesting  observation  is  that  the  ring  and  stringer 
stiffened  shells  already  have  a series  of  "reinforcing 
frames"  surrounding  the  cutout  consisting  of  the  rings  and 
stringers.  The  additional  stiffness  of  this  configuration 
requires  heavier  frames  in  order  to  alter  the  deformations 
significantly,  whereas  the  relatively  flexible  stringer 
stiffened  shell  does  not  require  heavy  frames  to  alter  the 
buckling  pattern.  In  order  to  give  a more  detailed  expla- 
nation of  this  effect,  consider  the  stiffening  effect  of 
the  ring  of  the  reinforcing  frame.  Since  the  volume  of  rein- 
forcing frame  remains  constant  for  a given  volume  ratio, 

W */W  , the  cross-sectional  area  of  the  frame  is  inversely 
sf  o’ 

proportional  to  the  distance,  d,  of  the  frame  from  the  cen- 
terline of  the  cutout  and  the  moment  of  inertia  to  the 
square  of  d.  Hence,  for  the  ring  and  stringer  stiffened 
shell,  the  added  stiffness  of  the  reinforcing  frame  is  large 
when  the  frame  is  along  the  cutout  edge;  but  for  the  frame  re- 
moved from  the  cutout  it  is  not  as  large  when  compared  to  the 
stiffness  of  the  existing  ring  stiffened  skin.  On  the  other 
hand,  for  the  stringer  stiffened  shell,  the  relative  increase 
in  stiffness  in  the  circumferential  direction  is  large  even 
for  a frame  removed  from  the  cutout;  since  without  the  frame, 
only  the  relatively  weak  skin  resists  bendinc  and  compressive 
loads. 


note  that  Figure  24  through  27  can  be  used  to 

the  required  combination  of  frame  volume  and  frame 

required  to  achieve  a certain  ratio  P /P  . 

cr  o 


de  termine 
location 


b 


i 
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For  example, 


consider  Figure  27  which  is  for  a ring  and  stringer  stiffen- 
ed shell  with  a 24"  x 24"  cutout.  Suppose  that  it  is  not 
possible  to  reinforce  the  cutout  at  the  edge,  but  can  be 
reinforced  at  position  2,  what  is  the  required  volume  of 
reinforcing  frame  to  achieve  a buckling  load  equivalent  to 
a similar  shell  without  a cutout?  The  required  volume  ratio 
of  the  reinforcing  frame  can  be  easily  determined  as  follows 
Draw  a line  through  P /P  = 1.0  that  intersects  the  frame 
position  2 curve.  A line  dropped  from  this  intersection 
to  the  axis  yields  the  required  volume  ratio  (and  therefore 
volume)  for  the  reinforcing  frame.  Note  that  this  proce- 
dure works  eaually  well  for  any  other  value  of  P /P  . 


IV.  CONCLUSIONS 

The  following  conclusions  can  be  made. 

(1)  A reinforcing  frame  equal  in  volume  to  the 
material  removed  for  a cutout  and  placed 
adjacent  to  that  cutout,  can  increase  the 
buckling  load  of  the  shell  to  or  above  that 
of  a shell  without  a cutout. 

(2)  For  either  a stringer  or  a ring  and  stringer 
stiffened  cylindrical  shell  with  rectangular 
cutouts,  the  optimum  position  for  a reinforc- 
ing frame  is  along  the  cutout  edges. 

(3)  Positioning  a frame  away  from  the  cutout  edge, 
in  most  cases,  drastically  reduces  the  buck- 
ling load  from  that  obtainable  by  placing  the 
same  volume  of  reinforcing  frame  adjacent  to 

' the  cutout. 

(4)  For  very  light  frames  there  may  be  a frame 
position  away  from  the  cutout  edge  which  is 
equally  effective  as  the  position  adjacent 
to  the  cutout. 

(5)  Frames  positioned  far  away  from  the  cutout 
ecce  msv  actually  reduce  the  oucklinc  load 
below  that  of  an  unreinforced  shell. 


A 


1 

(6)  It  is  possible,  with  the  use  of  STAGS  to  de- 
termine the  reinforcement  frame  cross-sec- 
tion necessary  to  efficiently  stiffen  a 
shell  with  cutouts  to  any  desired  strength. 
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Appendix  A 


Reinforcing  Frame  Geometric  Parameters 

The  geometric  Darameters  for  the  reinforcing  frames 
depend  upon  the  independent  parameters;  shell  stiffening, 
cutout  size,  volume  ratio,  and  frame  position  according 
to  Eqs.  (35-39).  For  completeness,  the  frame  distance, 
d,  cross-sectional  area,  Asf>  height  hsf;,  moment  of  inertia, 
Is^.,  and  torsional  stiffness,  GJsf,  of  the  reinforcing 
frame  are  listed  below  for  each  run. 
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Buckling  Loads  and  Buckling  Load  Ratios 

The  buckling  loads  and  the  ratios  of  P /P  and 

cr  o 

PC£./P » are  listed  below  versus  computer  run  number. 
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